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Abstract. Global existence for weak solutions to systems of nematic liquid 
crystals, with non-constant fluid density has been established in '27' and 12]. 
In this paper we extend the regularity and uniqueness results of Fanghua Lin 
and Chun Liu in 1221 for the systems of nematic liquid crystals (LCD). In I22| . 
the underlying system has constant density. In this paper, the regularity and 
uniqueness results are established for a density dependent LCD system. 



The flows of nematic liquid crystals can be treated as slow moving particles 
where the fluid velocity and the alignment of the particles influence each other. 
The hydrodynamic theory of liquid crystals was established by Ericksen [7] , [S] and 
Leslie [50], [H] in the 1960's. As Leslie points out in his 1968 paper : "liquid 
crystals are states of matter which are capable of flow, and in which the molecular 
arrangements give rise to a preferred direction" . In this paper we consider the 
simplified model for the flow of nematic liquid crystals: 



in Q X (0,T), where is a domain in i?", p : fl x [0,T] — >• K is the fluid density, 
p : fix [0,T] ^ R is the fluid pressure, m : x [0, T] is the fluid velocity and 

d : X [0, T] — > R" is the director field representing the alignment of the molecules, 
with n — 2,3. The force term Vd (g) Vd in the equation of the conservation of 
momentum denotes the 3x3 matrix whose ij-th entry is given by V^d • Vjd for 
1 < hj < 3. This force Vd(g)Vd is the stress tensor of the energy about the director 
field d, where the energy is given by: 



1. Introduction 



(1.1) 



+ V • (pu) = 0, 

{pu)t + V • (pu (g) u) + Vp = Alt - V • (Vd «) Vd) 
dt + u • Vd = Ad - /(d) 
V • M = 




and 



F(d) = — (|dp - l)^ /(d) = VF(d) = -(|d|2 - l)d. 
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In fact F(d) is the penahy term of the Ginzburg-Landau approximation of the orig- 
inal free energy of the director field with unit length. 

There is a vast literature on the hydrodynamic of the liquid crystal system. For 
background we list a few names, with no intention to be complete: [9], [13], |14) . 
[22], [23], [24], [4], [3], [4], [5], [6], [27] and [I^. Particularly, in [27] [12], the global 
weak existence of solutions to the flow of nematic liquid crystals was obtained for 
fluids with non constant density. In light of the regularity results to the pure fluid 
system established in [2 and [T9 it is natural to ask if the regularity results in [22 
can be extend to prove the regularity of the solutions for flows of nematic liquid 
crystals with non-constant fluid density. 

In this paper we focus on the regularity of solutions to the flow of nematic liquid 



crystals satisfying the initial conditions: 

(1.2) p{x,0) ^ po{x), < Ml < po{x) < A'h, 

(1.3) u{x,0) ^ uo{x), V • Wo = 0, Mo|an=0, 

(1.4) d{x,0) = do{x), \dn{x)\^l, 
and the boundary conditions: 

(1.5) u{-,t)\dn^O, d{-,t)\gn ^ do\gn. 



Existence of global weak solutions of (jl.ip . with the above specified data, has 
been established in [27] and [12] . In fact they have existence even without assuming 
the positive lower bound Mi. In [2^ and [T^, to derive global weak solution, a 
viscosity term eAp is added to regularize the first equation of the system (II. ip . This 
approach had been suggested in }25) . Our proof of regularity uses energy estimates 
introduced by Ladyzhenskaya on the approximate solutions. The added term in 
the second equation that results from the regularizing viscosity in the first equa- 
tion in |25j . contains the gradient of the fluid density. This term seems to create 
difficulties when it is used to establish the Ladyzhenskaya energy estimates for the 
approximate solution derived by the Galerkin method. Thus in the appendix we 
sketch a proof of existence for the global weak solutions to system (|l.ip without the 
introduction of the viscosity term for the density in the equation of the conservation 
of mass. In our case to obtain a classical solution we need to work with data that 
is more regular than the data used in [57^ and [T^ . 

We obtain interior regularity with a relatively weak conditions on the initial 
data. For more regular data, we are able to obtain solutions which are regular up 
to the boundary. In the rest of the introduction we briefiy describe our main results: 

Regularity in 2 dimensions: 

Theorem 1.1. Suppose that ^ is a smooth bounded domain in M?. Let pQ, uq and 
do satisfy fO)) - l775)) . Suppose that po e C^(f2), uq e H^{n) and do € H'^{Q.). 
Then, the system has a global classical solution (p, u, d), that is, Jar allT > 

and some a G (0, 1) 

we ci+"/2'2+"((o,r) X n) 
VpG c"/^'"((o,r) X n) 
de c^+"/^'2+"((o,r) X ii) 
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ci((o,r) X n). 

Corollary 1.2. Suppose {pQ,UQ,do) satisfies the conditions in Theorem 11.11 In ad- 
dition, Pa e Ci(0), Mo e C2+"(17) and do e C^+'^in). Then the solution {p,u,d)_is 
regular up to boundary, that is, the conclusions in Theorem 1 1 . 1 1 hold in[0,T] x $1. 

Indeed, more regularity on the initial data do implies more regularity on the 
boundary due to the second condition in (11.51) . Thus, yielding the regularity up to 
boundary, see [l5] . 

The regularity for the flow of the nematic liquid crystals in dimension 2 turns out 
to be not too difficult. We first establish the Ladyzhenskaya energy estimate (|2.19p 
and (|2.20l) . similar to that in [22]. Then we apply the regularity result for trans- 
port equations in [2 to obtain the Holder continuity of the fiuid density. Therefore 
Theorem 11.11 follows from the estimates and Holder estimates in [17] and a more 
or less standard bootstrapping between the three equations in the system ()1.1|) . 

Provided we have sufficiently small data or we work with sufficiently short time 
in 3D, we also obtain interior regularity. Given more restriction on data, the regu- 
larity can be obtained up to the boundary. 

Regularity in 3 dimension with small data or short time: 

Theorem 1.3. Suppose that ft is a smooth bounded domain in R"^. Let pQ, uq and 
do satisfy U.2\} - il.5\} . Assume that po € C^(ri), uq e H^{fl) and do £ H'^{fl). 
Suppose that {p,u,d) is a weak solution to the system (jl.ip in Theorem \5.1\ Then 

1. There is a positive small number eo such that if 

(1-6) \\uo\\Hi(n) + l|Vcfo|lHi(o) < eo, 

then the system lll.l]) has a classical solution (p,u,d) in the time period (0,T), for 
all T > 0. That is, for some a G (0, 1) 

u e ci+"/2'2+"((o,r) X n) 
Vpe c"/2'"((o,r) X n) 

(1.7) , 

ci+"/2^2+"((o,r) X fi) 
peC\iO,T) X n). 

2. For general data, there exists a positive number 5o = So{po,uo,do) such that 
|j.7[ ) holds in the interval (0,T) for T < So- 

Corollary 1.4. Suppose in addition to the hypothesis in Theorem 11.31 that po G 
C^{n), uo € C2+"(17) and do e C2+"(f]), then the solution is regular up to the 
boundary for data small in the sense (|1.6p or for large data and sufficiently short 
time. 

The proof of the regularity of the solution to the system (|l.ip in dimension 3 
takes the same approach as in dimension 2 but is much more complicated. First 
in contrast to the cases of dimension 2, we only get the Ladyzhenskaya energy 
estimates when either the initial data is small in the sense as described in (|1.6|) 
or T is small. Our calculations and estimates are based on those in [22], with 
interesting modifications. We use ideas of [22] making it work in a rather different 



4 



MIMI DAI, JIE QING, AND MARIA SCHONBEK 



way. We keep the potentially small terms ||u||i2 and ||V(i||L2 instead of throwing 
them away. This gives a more unified way to derive the Ladyzhenskaya energy 
estimates in the cases: 

• of small data 

• for short time. 

After having the Ladyzhenskaya energy estimates, in contrast to the two dimen- 
sional case, we do not have the Holder continuity for the fluid density. Instead we 
observe that we have small oscillations of the density over small balls in 51 x [0,r] 
provided that either the initial data is small or for short time. This turns out to be 
enough to carry out the frozen coefficient method to improve the regularity of the 
fluid velocity. We refer the reader to [TH] for a reference of the frozen coefficient 
method. We give the idea of the method in Appendix 6. Our key lemma on the 
oscillation of the fluid density is as follows: 

Lemma 1.5. Suppose that Q is a smooth bounded domain in . Let pQ, and 
do satisfy U.2]} - U.5]} . Assume that po £ C^{Cl), uq G H^{il) and do e H'^{il). 
Suppose that [p, u, d) is a weak solution to the system (jl.ip in Theorem I5.il Let 
ti £ (0,T) and p G Vt, define 

Ap.ti) = {BrM n r!) X ([ti - ro,ti +ro] n [0,T]). 

Then, for any e > 0, there exists eq > and tq > such that for p £ fl and all 
T > 0, 

(1.8) sup \p{q,t2) - pip,ti)\ < e, 

(ij,t2)6-4(p,tj) 

provided that either 

ll"o||l/i(j2) + ||Vdo||H2(o) < eo orT < do- 
Remark 1.6. The interior regularity in Theorem 11.11 and Theorem 11.31 is obtained 
by bootstrapping argument. 

To close this section, we state the uniqueness of solution in the following sense. 

Theorem 1.7. Let {p,u,d) he the solution to system iLl]) and iL2\) - !!775\) obtained 
in Corollary \L2\ for two dimensions or in Corollary \1.4\ for three dimensions. Let 
{p,u,d) be a weak solution to system with iL2\) - n~5\) satisfying the following 

energy inequalities: 

(1.9) / \p\^dx< f \po\^dx, 
Jn Jn 

(1.10) [ Ip\u\^ + ^\\7d\^ + F{d)dx+ f f \\7u\^ + \Ad- f{d)\^dxdt 
JO 2 2 Jo Jfi 

- X + ll^M^ + F{do)dx. 

Then (p, u, d) = (p, u, d) . 

Remark 1.8. Note that the weak solutions established in [27] and [12] satisfy the 
inequalities recorded above. Hence these solutions coincide with solutions we ob- 
tained in Corollary 11.41 due to the uniqueness Theorem 1 1.71 
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2. Classic solutions to Nematic Liquid Crystals Systems in 

In this section we use the Ladyzhenskaya energy method |18j to show that u G 
i-([0,r],i/i(f^))ni'([0,T],i?2(f^))anddGi-([0,r],ij2(C!))nL2([0,T],i?3(f^)) 
Then we apply a resuh from 2 to get the Holder continuity of the density p. With 
the continuity of the density p we apply the so-called frozen coefficient method to 
get higher regularities for p,p,u,d by bootstrapping. The proof of the Ladyzhen- 
skaya energy inequality is similar to that in [22' . The key inequalities used often in 
this paper are the following Gagliardo-Nirenberg inequality (cf. 10 ) : 

Lemma 2.1. (Gagliardo-Nirenberg) If fl is a smooth bounded domain in M", then 

(2.11) \H\hin) < t^lkllL(o)(l|V^||i2(o) + Mlh^n)), 
when n — 2 and 

(2.12) Ml.^n) < C^hllL^(o)(l|Vzi||i.(f,) + h||i2(f,))i, 
when n — 3. Moreover, = 0, then 

(2.13) 1^1114(0) <C||w|li.(o) II V^illi.(o), 
when n = 2 and 

(2.14) <C|klU^(0)l|V^i||i.(f,), 
when n = 3. 

2.1. The L°°{H^) and L'^{H^) estimates of the velocity. Our strategy is the 
same as in j22j . We first establish the desired bounds for the Galerkin approximating 
solutions {pm,Um,dm) in the sequence that one has from Galerkin method when 
proving the existence (cf . Section 5 of the appendix in this paper) . By passing to 
the weak limit we then obtain the desired bounds for the weak solution {p,u,d). 
First we set 

(2.15) 'i>l{t) = \\Vun\h + \\Adn\h. 

Then to show that m™ G L°°{H^)f]L^{H^) and rf" G {H^) f] {H^) wc cal- 
culate -^^m- Namely, 

2 dt Jq Jq 

= - [ [ |VAd™|2dx 

Jn Jn 

+ [ VAd™ . V(ii™ . Vd™) - Ad" . A(/(d"))dx 
Jn 

+ [ -p™(u™ . Vu™)mJ" - wT'Vd^Ad^dx 
Jn 

= - [ /9">J"|2da;- [ \VAd"'\^dx 

Jn Jn 
+ I + II + II + IV, 

here we used the equation (|5.53l) . We find 

/ < e||VAd"||2, + C\\Wu"'\\l4Wd"'\\l, + C|lu™|li4|lAd'"||i4 + Cllu™!]^^ 
< e\\Au"'\\l2 + 2e||VAd™||i2 + + + + 
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by Gagliardo-Nirenberg inequality, the basic energy estimate (|5.59p and the fact 
that 

llvV"|U4 <c(|lAd'"|U4 + |id-|U4). 

And similarly, 

/// < eWlSu^Wl. + e\\uT\\l. + C< + C$^, 

IV < e||VAd'"||i. + e\\uT\\h + C^tn + C^'m- 
Moreover, from equation (|5.53p . we have 

(2.16) II A«"||i. < 2M2\\uT\\h + ZllVArf^lli. + C< + + C<i>l 
and 

(2.17) Mill^^rili^ < ZllAu-lli. + 2||VAd'"||i. + C< + C*^ + C*^. 
Therefore, we arrive at 

(2.18) j^<i>lM + \\Aun\h + llVAd-lli. < C<(t) + C, 

which implies, as in [22; , again in the light of the basic energy inequality (j5.59p . 

^l{t) < (<f>'(0) + C7)eC^/o"(C(*)+i)'i* < ($(0)2 + C)e^T+c 
where ^^(o) = ||Vuo||i2 + ||A(io|||2, and 

/ ||Au"||2,dt+ / ||VAd"||2,df < C. 

Hence we have obtained 

Theorem 2.2. Suppose that fl is a smooth bounded domain in and that $(0)^ = 
||Vuo||^2 + IIAdoll^j < oo. Suppose the initial data {po,uo,do) is as in The- 
orem Then there exists a constant C such that the approximating solutions 
(p*", u™, d"*) of the system (|1.1[) obtained in Theorem \5.1\ satisfies 

(2.19) llVu^ll^, + ||Ad™||2, < e^^+^($2(o) + C) 
for all t G [0, T] and 

(2.20) [ \\Au"'\\l2dt+ [ ||VAd™|||2dt < C. 
Jo Jo 

From the last theorem it follows that 

Corollary 2.3. Under the same hypothesis of last theorem, there exists a solution 
{p,u,d) to (|l.ip which satisfies the energy inequalities 

(2.21) l|Vu||i2 + ||Ad||2, < e^^+^($2(0) + C) 
for all te[0,T] and 

(2.22) / \\Au\\l2dt+ [ \\\/Ad\\l2dt<C. 
Jo Jo 

Proof: It follows by extracting a subsequence of the Galerkin approximations 
(p™, u™, d™) and passing to the limit. 
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2.2. Holder continuity of the fluid density. Next we recall a regularity lemma 
for the transport equation from 2] to get Holder continuity for the fluid density p. 

Lemma 2.4. ([2J Suppose that u e L°°{[Q,T],H^{VL))f^L^{[Q,T\,H^{Vt)). And 
suppose that 

pt + u ■ Vp = 

in rt X (0,T). Assume that p{0) G C^{fl) and that C is smooth and bounded. 
Then p € C^ifl x [0, T]) for some a G (0, 1) which depends only on the initial data, 
T and f2. 



2.3. Proof of Theorem 11.11 The proof is more or less standard, particularly 
after the work in 2 . We sketch a proof here for completeness. Rewrite the third 
equation in (jl.ip as 

(2.23) dt- Ad^ -u-Vd- f{d). 

It yields from the basic energy inequality (|5.51|) and estimate (|2.19|) that , 

u-Vde L°°{0,T;L''{n)), for any r e [l,oo). 

Then due to standard estimates of solutions to parabolic equations (see [T^ and 
[I]), from the equation (|2.23p we have 

deW^'^W^''), for any r,s e [l,oo), 

which implies that 

(2.24) de c"/2'i+"([o,r] X f7) 

for any a £ (0, 1). Now we go to the second equation in (jl.ip and use the above 
estimates for d to improve the estimate on u via the frozen coefficient method |19) , 
after we have the Holder continuity for the fluid density p, as done in [2] to derive 
that 

for some q > 1 and any p > 1. Therefore we have 

u e c"/2^"([o,r] X n). 

Next going back to the third equation in 11.11 by the standard Holder estimates for 
parabolic equations (see [T7] Theorem 5.1 in Chapter VII) we have 

d e C^+'i'^+"{{0,T) X n), 

for some a G (0, 1). Therefore we are able to go back to the the second equation in 
()1.1|) . again via frozen coefflcient method [19], to derive that 

u g Ci+t^2+"((o,T) X fl). 

Finally p E C^((0, T) x ft) follows from the regularity of u and the regularity of the 
pressure p follows easily from the regularity of (p, u, d) similarly as in [2]. 

□ 

Remark 2.5. By bootstrapping argument, we can get higher regularity. 

Proof of Corollarv 11.21 follows by Krylov's Theorem 10.3.3 in [15 , which for 
convenience of reader, we recall here [15] . 
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Theorem 2.6. For any h e C"/2'"((0,r)xfi) and boundary g e C^+°'/^-^+"{{0,T)x 
n) there exists a unique function u G C"'^+"/^'^+"([0, T] x 57) satisfying the equation 
Lu — Ut — h in (0, T) x 57) and equal g on boundary of (0, T) x 57). 

We apply this to both the second and third equation in (jl.ip . Applying to the 
second equation, h = u ■ Vd + f{d), g — do and L = A. Applying to the third 
equation, h = u - Vw + Vp + V • (V(i(Xi Vd), g — Q and L = A. Proceeding as in the 
last theorem yields h e C"/^'"((0,T) x 57) for both equations. On the boundary 
of 57, u = 0, d = do and do G C'^+"(57) by hypothesis. Hence the conclusion of the 
Corollary follows. 



3. Classical Solution to Nematic Liquid Crystals System in 3D 

In this section we consider the solutions of (jl.ip for 57 e R'^ a bounded domain. 
We establish regularity in two cases 

• global regularity with small initial data 

• short time regularity. 

First we adopt the idea from [22] to derive the Ladyzhenskaya energy estimates. 
In contrast to the cases of dimension 2, we will only get the Ladyzhenskaya energy 
estimates for the above two cases as in [52]. Our calculations and estimates are 
similar to those in [22j . with interesting modifications. After having the Ladyzhen- 
skaya energy estimates, unlike the case of dimension 2, we will not have Holder 
continuity for the fluid density right away. Instead, we observe that we can have 
the oscillation of the density over small balls in 57 x [0,r] to be small, provided 
that either the initial data is small or we work for short time, which turns out to 
be enough to carry out the frozen coefficient method to improve the regularity of 
the fluid velocity. 



3.1. Ladyzhenskaya Energy Estimates. Our derivation of the Ladyzhenskaya 
energy estimates in p2] in dimension 3 is rather an interesting modification of 
the original Ladyzhenskaya's argument for the pure fluid systems. In [22 for the 
argument to work, it needs either the viscosity for the fluid to be very large or the 
time to be very short. We use the same idea, but, instead we assume the initial 
data to be small or the time to be short while the viscosity of the fluid is a fixed 
constant. For the convenience of the arguments in our context, without loss of 
generality we take the constant to be 1. Let us set as before: 

(3.25) <l>2„(t) = ||Vu™||i. + ||Ad"||i. 

Then, as in the previous section, we will first derive the Ladyzhenskaya energy 
estimates for the Galerkin approximate solutions (p™, u™, d™) and pass to the weak 
limit to obtain the Ladyzgenskaya energy estimates for the weak solutions (p, u, d) 
to the system as desired. 
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Again, Using u"^\ga = d"'-\ga = 0, and Ad"^\ga = 0, by integration by parts, it 
still follows that, 

~^l^(t) = / Vii" • \7ul"dx + [ Ad™ • AdV'dx 

! Jn 
VAd™ • V(ii™ • Vd") - Ad" • A{f{d"'))dx 

-p'^iu"' ■ Vii™)w7' - u'l'\7d"'Ad"'dx 

= - [ / |VAd™|2dx 

+ I + II + II + IV, 

We will proceed the same way as we did in dimension 2 except the Gagliardo- 
Nirenberg inequality is different in dimension 3 from that in dimension 2. More 
importantly we will keep the terms and ||Vc?™||l2 whenever necessary. We 

may derive 

/ < e\\Au^\\l. + e||VAd"||i. + + || Vd™||i.)(C< + C), 

//<C$2„ + q|Vd"|U2, 
/// < e\\uT\\h+e\\Au"^\\h + C\\u"Tl^^I + C-fl, 

and 

IV < £||u™||2, + £||VAd"||2, + C||Vd"||i2$^ + C||Vd™||2,. 
To relate Am™ back to u™, from equation (I5.53p . we have that 
II AM™||i. < 2Ah\\uT\\h + 2||VAd™||i. + (||u™||i. + ||Vd™||i.)(C^< +C) + C^l 
and 

M,\\uT\\h < 2|| Au™||i. + 2||VAd™||i. + (||M™||i. + ||Vd™|li.)(C<I>f„ +C) + C^l 
Therefore we arrive at 

II Au™||i. + ||VAd™||i. + j^^l < C<^1 + (||u™|U. + ||Vd™|U.)(C^*^. + C). 

Set 

^™-*™ + lh^"llL^ + l|Vd'"||L2 

and observe that 

(3.26) j^^l < (C + C(||«™|U. + II Vd™|U.)|.^)|.^ 
or set 

*™ = <I'™ + ll""llL^ + llVd™|U. + l 

we have 

(3.27) < 

It is clear from p.27p one can prove the Ladyzhenskaya energy estimates when T 
is small. To get the Ladyzhenskaya energy estimates for small initial data we will 
use an idea similar to that in [22]. Suppose that 

hofm + WdoWl^^eo. 
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Recall that by the basic energy estimate we have 

\\u^\\l2 + ||Vd™|U2 < \\uo\\l^ + WdaWl, < 00. 
Then we claim that, if is so small that 

(3.28) C9o{ie'^+^9of < 1, 
then 

for all t. First we prove the claim for t G [0, 1]. Assume otherwise, there must be 
to G (Oil) such that 

(3.29) ^lito) ^ Ae'^+'eo 
and 

for all t G (0,to]- Therefore, from p.26p . by the choice of 9o in p.28p . we have 

il$2 < (c + 1)1,2 

for all t G (0,^o) and ^^^(0) < 29o, which implies that 

^lito) < e^+i^^(0) < 2e^+% 

and thus contradicts p. 291) . For i > 1, we simply observe, as in [22 , that the basic 
energy inequality (|5.5ip 

^lit)dt< 200 

L 

implies that there is to G (^ ~ 1 1 1) such that 

Then one may repeat the above argument to conclude that 

As in the case of 2 dimension, passing to the limit the existence of weak solutions will 
follow from the uniform estimates we have obtained. This weak solution satisfies 

Theorem 3.1. Suppose that il is a smooth bounded domain in M.^ . Let pQ, ug and 
do satisfy lll.S\) - l[T75\) . Assume that 

\\uo\\h^ + WdoWn^ < oo. 

Let {p,u,d) be a weak solution to the system (jl.ip with data {po,Uo,do). There is 
Co > such that if 

Ikolki + II^oIIh^ < eo, 

then 

(3.30) \\Vu\\l, + \\Ad\\l,+ [ {\\Au\\l. + \\S/Ad\\l.)dt<Ci\\uo\\Hi + \\do\\H-) 

Jo 

holds for any positive T. For data with no smallness condition there is a 6o de- 
pending on the data such that \3.3U\) holds for all T < So- 
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3.2. Oscillation of the Fluid Density. In this subsection, we show that the 
oscillation of the density p over small balls in space time. This estimate follows 
by the estimate on u from the Ladyzhenskaya energy estimates obtained in the 
previous subsection. For convenience of the readers we restate Lemma 11.51 



Lemma 3.2. Suppose that Q is a smooth bounded domain in K.^ . Let pa, uq and 
do satisfy U.2 \) -S1.5 \) . Assume that po e C^{0,), uq £ H^{n) and do G H^{n). 
Suppose that {p, u, d) is a weak solution to the system (jl.ip constructed in Theorem 
15. Jl ( see Appendix). Let ti S (0, T) and p G Q, define 

-4(p,ti) = (Broip) nn)x {[ti-ro,ti+ ro] n [0, T]). 

Then, for any e > 0, there exists eo > and tq > such that for p £ and all 
T > 0, 

sup \p{q,t2) - p{p,ti)\ < 
(9,t2)e.A(p_tj) 

provided that either 

Ikollffi(o) + l|Vdo|lH2(a) < eo or T < So- 
Proof: The equation of the conservation of mass is 

(3.31) p^ + u-Vp = 0. 

For reasonably regular velocity we may solve p.3ip by the method of characteristics. 
Due to Osgood theorem [29], since u S for almost all t S [0,T] and in light of 
the Ladyzhenskaya energy estimates in Theorem 13.11 there is a unique solution to 
the Cauchy problem corresponding to (I3.3ip by finding trajectories of the liquid 
particles 

^=u(2/,t), y\^^t=x- xen, te{0,T) 

and defining 

(3.32) p{x,t)^ po{yiO,x,t)). 

Step 1 Fix a time t G [0,1"]. Let xi and X2 be two arbitrary points from 
satisfying \xi — a;2| < d < 1. For any r € (0,t), assume 

yi = y{T,xi,t), y2 = y{T,X2,t), 

then the difference z{t) = yi — y2 is the solution to the Cauchy problem 

dz 

— = u(yi,T) - u(?/2,t), z\r=t^Xi-X2. 

dr 

For < a < 1/2, by a standard Sobolev embedding theorem in |TU], we have 

^<C\\u{T)\\H^n)\zr, re(0,t). 
Integrating form r to i, 

(3.33) |z|i-" < \xi-X2\'-''+CT-^\\u\\mo,T-HHn)). 

As a consequence of Theorem 13. 11 we know that \z\ is as small as needed provided 
that \xi — X2\ and ||wo||/fi + Udollff^ are small or T is small. Therefore 

\p{xi,t) - p{x2,t)\ = \po{y{0,xi,t)) - po{y{0,x2,t))\ < C\z\ 

is as small as one wants provided that \xi — X2\ and ||uo||hi + IIc^oIIh^ are small or 
T is small. 
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Step 2 Fix a point x £ H,. Let ti,t2 E [0,T] arbitrary, let 

yi = y{T,x,ti), y2^y(T,x,t2). 

Assume that x' — ?;2|r=ti- Then due to uniqueness, the integral curve y2{'T) can be 
considered as a solution to the Cauchy problem 

■"(2/2, t), ?/2|r=ti = x' 



dy2 



dr 

with initial data at r = ii. Hence, the difference z(t) — yi ~ y2 is the solution of 
the Cauchy problem 

dz 

— = u{yi,T) - u{y2,T), z\r=ti^x-x'. 
dr 

By the definition of a;', we have 

u{y2{s,x,t2),s)ds. 
Therefore, again due to a standard Sobolev embedding Theorem, 

|x-a;'|< / \u{y,s)\ds<C \\u\\H'-'{n)ds < C\\u\\L2(^Q^T-M^n})\ti - t2\"^ ■ 
Jti Jti 

By Step 1, we conclude that 

\p{x,h)-p{x,t2)\ = \po{y{0,x,ti)ypo{y{0,x,t2))\ = \po{y{0,x,h))-po{y{0,x' ,h))\ 

is as small as needed provided that \ti — t2\ and ||7io||/fi + ||(io||/f2 are small or T is 
small. This completes the proof of the Lemma. 

□ 



3.3. Proof of Theorem 11.31 In this subsection, we finish the proof of Theorem 
11.31 This is done by combining Theorem 13.11 Lemma 11.51 from the previous subsec- 
tions, the frozen coefficient technique applied to LP{L'^) estimates, Holder estimates 
and a bootstrapping argument between the three equations of the system (jl.ip . 
In the appendix wc briefiy describe the frozen coefficient method. For simplicity, 
we show how the method works for the density dependent NSE and obtain the 
estimate for L'^ space with q > \/3. We apply now this method to our approxi- 
mating solutions (p"*, M™, d™). These solutions for each m satisfy the conclusion 
of Theorem 11.31 We now show that the estimates are uniformly in m applying the 
frozen coefficient method, in which case we let q = 2. That allows us to pass to the 
limit (p, u, d) and obtain (jl.7p for weak solutions. We now show briefiy the steps 
to yield the necessary uniform bound. 

Here to simplify the notation we denote the approximating solutions by {p,u,d). 
We notice first that since 

ue L°°(0,T;L^(fi)), Vde L°°(0,T;L6(f])), 

we have 

u-Vde L°°{0,T;L^in)). 

By standard parabolic estimates on the third equation in (jl.ip (cf. [TT^ and [Tj), 
we have 

deW^'PiW^'^), 
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for all 1 < p < cx), which implies that Vd £ L°°{0,T; L'^{i})) for any q £ (l,cxo). 
Thus, we have 

li- Vde L°°(0,T;L9(f7)), Vge(l,6). 
Applying the same standard parabolic estimate on the third equation in (jl.ip yields 

(3.34) d e VF^'f (T/F2'«), Vp e (1, oo) and q e (1, 6), 
which implies that d £ C"/2,i+"([o, T] x 0.) for some a e (0, 1) and 

WdAde L°°{0,T;L''{n)), Vqe(l,6). 

In the second equation of (jl.ip , the estimates for the conservation of momentum 
with constant density can be extended to the non-constant density cases when 
Lemma 11.51 is available. This is done via the frozen coefficient method. 

We know that 

u-Wue L°^{0,T;Li{n)). 
Now we apply the frozen coefficient method using the oscillation estimates for the 
density, fLemma ll.51) to yield 

Thus u e L°^{0,T;W^'^{n)) and u • Vu e L°°{0,T; L^{n)). Repeating the above 
argument yields 

(3.35) uGW^'PiW"^'^), Vpe(l,oo), 

from where it follows that u G C"(rj x [0,T]). 

Back to the third equation in (jl.ip . we conclude that 

c'i+f'2+"((o,T) X n) 

for some a G (0,1) via the standard Holder estimates for parabolic equations (cf. 
[17) . Chapter VII). From here the argument for the pure fluid systems works with 
no further significant modifications. This completes the proof of Theorem 11.31 

Proof of Corollarv 11.41 Follows by Krylov's Theorem [16] just as in the case of 
two dimension. 

4. Uniqueness of Solution 

In this section we establish Theorem 11.71 For the LCD system with constant 
density, Lin and Liu [22] proved that the solution (u, d) is unique provided u G 
L°°(0, T; iJi) and d S L°°(0, T; H'^). The idea is to calculate the energy law satisfied 
by the difference of two solutions and establish a Gronwall's inequality. In our case, 
to calculate the energy law of the difference of two solutions it has some extra terms 
involving the density. Hence the estimates are more involved requiring additional 
bounds on the strong solution {p,u,d) to yield a Gronwall's inequality. In 2D, we 
need 

(4.36) Vp,Vue L°°{{0,T)xn), ut,u -Vu £ L°°{0,T;L'>{n)),q > 2. 
In 3D, we need 

(4.37) V/9,Vu e i°°((0,r) X rj), ut,u-Vue L°°{0,T;L^{n)). 

With the assumption on data, po G C^i^l), uq G C^+^'iO.) and do G C^+^iCl), 
the solution {p,u,d) from Corollarv 11.21 or Corollarv 11.41 satisfies (|4.36p or (|4.37p . 
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respectively . 

Proof of Theorem 11.71 First recall that the solution {p,u,d) from Theorem 11.11 
and Theorem 11.31 satisfies energy equality: 

(4.38) [ Ip\u\^ + ^\\/d\^ +F{d)dx+ [ [ \\/u\^ + \Ad- f{d)\^dxdt 

= \po\uo\'' + ^iVdol' + F{dQ)dx. 

The density p is the strong solution of the transport equation, hence it satisfies that 

/ p^dx — I p^dx. 
Jn Jn 

On the other hand side, p is a weak solution of the transport equation and Mi < 
p < We have by hypothesis that 

(4.39) / p^dx < / pldx. 

Jn Jn 

Thus, 

(4.40) i / \p-p\^dx^^ [ p^dx+]- [ p^dx- I ppdx 
^ Jn ^ Jn ^ Jn Jn 

< / p^dx — / ppdx. 
Jn Jn 

Since p S C^{[0,T] x Cl), we can take p as a test function. Thus, multiplying 

Pt + u ■ Vp — 
by p and integrating by parts yields 

(4.41) / p^dx — / ppdx ~ — / pptdsdx —I {u ■ S/p)pdsdx 
Jn Jn Jo Jn Ja Jn 



I {u ■ \/p)pdsdx — / {u ■ \7p)pdsdx. 
Jo Jn Jo Jn 



lo Jn Jo Jn 

Here we used again that p is a classical solution of the transport equation. Substi- 
tuting the equality (|4.41l) in (j4.40p gives 



(4.42) 



- / \p — p\^dx< j j p{u — u)Vpdsdx 
2 Jn Jo Jn 

(p — p){u — u)W pdsdx. 



Jn 



Next, calculate the following term 

1 ( p\u~u\^dx + ]- [ iVd-VJpdx 

2 Jn 2 Jq 



n 

- / {p — p)\u\'^ dx + — f p\u\'^dx + — I p\u\'''dx — j pu-udx 



1 ( \\/d\^dx + l- [ iVdpdx- [ \7d®\7ddx, 

2 Jo Jn Jn 



where Vd Vc? denotes the 3x3 matrix whose ij-th entry is given by V^d • Vjd 
for l<i,j <3. 



REGULARITY OF SOLUTIONS TO THE LIQUID CRYSTALS SYSTEMS IN AND 15 



Using energy equality (|4.38p for the regular solution {p, u, d) and inequality ()1.9 
for the weak solution (p, u, d) combined with the last equality yields 

(4.43) i [ p\u~u\^dx+l- I iVd-Vdpdx 

< - / / \^u~Vu\'^dxdt- / / |Ad- AJpda;dt 
Jo Jn Jo Jn 

+ ^ / {p - p)\u\'^dx ~ I pu(^u- po\uo\'^dx ~ 1 Vrf(g) Vd" iVdopda; 



2 



t 



F{d)dx- [ F{d)dx + 2 [ F{do)dx- [ [ \f{d)\^dxdt 
Jn Jn Jo Jn 

[ [ \ f{d)\^dxdt-2 [ [ Vu-Vudxdt-2 [ [ Ad ■ Addxdt 
Jo Jn Jo Jn Jq Jn 

'iff Ad - f{d)dxdt + 2 f f Ad-f{d)dxdt. 
Jo Jn Jo Jn 



Since u,d & C"'^+"/^^^''""([0, T] x Q), we can take u, d as test functions for the weak 
solution u, d. Thus it follows that 



(4.44) / pu u — pq\uq\'^ dx — / / puutdxdt + I I pu[u ■ S/u)dxdt 
Jn Jo Jn Jo Jn 

/ pu{u ■ \7u)dxdt — / pu{u ■ \i'u)dxdt, 
>o Jn Jo Jn 



(4.45) / Vd<»Vd-\Vdofdx =^ - f f Addtdxdt + f f Ad{u-Vd)dxdt 

Jn Jo Jn Jo Jn 

- f f Ad- Addxdt + f f f{d) Addxdt. 

Jo Jn Jo Jn 



Indeed, formally, to get (I4.44p . we multiply equation 

{pu)t + div{pu (g) u) + Vp = Am - V • (V(i VJ) 
by u and integrate by parts. To get (|4.45p . we multiply equation 

dt+u-\/d = Ad~ f{d) 
by Ad and integrate by parts. 
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Substitute (|444l) and (|445l) in (|443l) . Add (iT42|) combined with the equahties 
in equations yields 

(4.46) 



\ I \P- P?dx + i I p\u- u\^dx + \ j \^d- \7d\^dx 

< - /" /" |Vu- Vwl^dxdi- [ [ \Ad- Ad\^dxdt 
Jo Jn Jo Jn 

+ / {p — p){u — u)\i' pdxdt — / / pS/ u\u — u\'^ dx 
Jo Jn Jo Jn 



[p- p){u-u){ut+u-Vu)dxdt+ I I u{Vd-'Vd){Ad- Ad)dxdt 
Jn Jo Jn 

[ {u-u){Vd^Vd)Addxdt+ [ [ {Ad- Ad){f{d) - f{d))dxdt. 
Jn Jo Jn 

Recall that the regular solution {p,u,d) satisfies (14.36^ in 2D and (|4.37l) in 3D. 

By Holder and Gagliardo-Nirenberg inequalities on the right hand side of ()4.46p . it 

follows that 

(4.47) i ^ \p{t) ~ mi' + mwit) - uit)\' + ivd(<) - vmi'dx 

<C [ [ \p- pf + \u-u\^ + \Wd-Wdfdxdt. 
Jo Jn 

To handle the last integral in (|4.46p . we used the fact that < 1 and hence 

\f{d) - f{d)\ < C\d - d\ by the definition of f{d). Thus, 

/ / \f{d)- f{d)\^dxdt<C [ [ \d-d\^dxdt<C{n) [ [ \Vd-Vd\^dxdt 
Jo Jn Jo Jn Jq Jn 

where the constant C depends on space domain not on time T, and C depends 

on the dimension of the space. 



Using the lower bound of p > Mi > 0, and Gronwall's inequality to (I4.47|) yields 
i \p{t) - Pit)? + pWit) - u{t)\' + \Vd{t) ~ Vd{t)\^dx 

< [ |p(0)-p(0)p+p(0)|w(0)-u(0)p + |Vd(0)-VJ(0)pda;e^* 

Jn 
= 

for alH > which implies 

p — p = u — u = d — d=0. 
This completes the proof of Theorem 11.71 



5. Appendix: Existence of Weak Solutions 

In this appendix we sketch an existence theorem for Galerkin approximations. 
As mentioned in the introduction the existence of global weak solutions to the flow 
of nematic liquid crystals have been established in [37] and in [T2] for non constant 
fluid density. Unfortunately the Ladyzhenskaya energy estimates do not seem to 
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work for the Galerkin approximate solutions constructed in |27j and in [12j . The 
Galerkin approximate solutions constructed here will possess the Ladyzhenskaya 
energy estimates. When the initial fluid density has a positive lower bound, we 
are able to derive estimates on and so that we can employ the compactness 
lemma of Lions- Aubin. Since the Galerkin method has been widely used for fluid 
systems as well as on the system of liquid crystals we will be rather brief (cf. [2] , 
[22], [22], m)- 



5.1. Galerkin Approximate Solutions. We construct a sequence of Galerkin 
approximating solutions that are uniformly bounded. These are the solutions used 
in Section 2 and Section 3. The bounds obtained there through Ladyzhenskaya 
method yield a subsequence that will converge to the classical solution. 
Let us first state the existence theorem for global weak solutions: 

Theorem 5.1. Assume that uq £ L^(fi) and do G H^{Q) with dolan & H^/'^{n). 
The system lil.l}) with initial boundary conditions il.S\) - !!775\) has a weak solution 
(p™, u™, d™) for each m — 1,2, 3, ... satisfying, for all T G (0, oo), 

(5.48) < Ml < < M2 

(5.49) It" e L^{0, T; H^{n)) D L°°{0, T; L^{Vl)) 

(5.50) d'"" £L^{0,T■H^{n))C^L°°{Q,T■,H^{n)), < 1, 

and the energy inequality 
(5.51) 

/ (|Ad'" - /(d")P + \Vu"'\^)dx + ^ / (JpIu™!' + JlVd^p + F{d"'))dx < 0. 
Jn dt Jq 2 2 

Proof. Let 

nin) = closure of {/ e C^{n,W^) : V • / = 0} in L^{n,R^) 
and {4>i}°^i be an orthonormal basis of TL and satisfying: 

-A(j), + \/P, = in n 
0j = on dn 

for i — 1,2,.... Here < Ai < A2... < A„ < ... with A; — > 00. In other words, 
we choose an orthonormal basis of H which consists of the eigenfunctions of Stokes 
operator on CI with vanishing Dirichlet boundary condition (see [31 ). Let 

P,n ■ U Hm = span{(f)i, ...0„J 

be the orthonormal projection. We seek approximate solutions {p™,u™,d™) with 
e Tim, satisfy the following equations: 

(5.52) + • Vp'" - 0, p'"(x,0)=po(x) 

(5.53) Pm(p™^u") = Prn(Au™ - p™u" • Vu™ - V • (Vd" (g) Vd"M) 

ot 

(5.54) d't" + • Vd" = Ad" - /(d'") 

(5.55) M™(x,0) = Pmuoix), d™(x,0) = doix), d'^ix,t)\an = do(a;)|ao. 
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(5.58) 



Let 

m 

(5.56) u^ix,t)^Y.9T{t)M^), 

1=1 

with 5™ (t) e [0, T] . Hence (|5.53p is equivalent to the foUowing system of ordinary 
differential equations: 

m , m m 

(5.57) J2-^T'{t)-9r{t) --J2^Tk'it)9r{t)9Tit) ~Y.Cigr{t)+V"-^^{t), 

2—1 i,k i—1 

for j — 1, 2, m, where 

Here <pj{x) is the kih component of the vector (pj{x). And 

rn „ 

w'"(-,0) = V5r(O)0»(a^), where g^iO) = uoix)Mx)dx. 
^=l 

To finish the proof, we need the following lemma. 

Lemma 5.2. There exists a weak solution (p^^ , u"^ , d"^) to the problem \5.52]) - 
\5. 55\) in Qt = 57 x [0,r], for any T e (0,oo), satisfying 

Ml < /o" < Ma 

e L°°{0, T; L^) n L2(0, T; H^) 

d"' £ L°°{0,T;H^)r]L^{0,T;H^), \d"'\ < 1. 

Moreover, (p™, u™, d™) is smooth in the interior of Qt and satisfies the energy 

equality, 

(5.59) 

/ |Vw'"p + |Ad™-/(d™)|2da; + -^ [ (i(p'"|u'"|2 + |Vrf™n + F(d"))da; = 0. 
Jn dt Jq 2 

The proof of this lemma is based on an application of the Leray-Schauder fixed 
point theorem. Let = T,T=i^T'^^ ^ C\Q,T;'H,n)- For each m let be a 
solution to 

(5.60) p + w"-Vp = 

with initial condition p(-,0) — pq. Let c?™ be a solution to 

dt + • Vd = Ad - /(d) 

with initial condition d(-,0) — dp and boundary condition d{x,t) = do{x). The 
reason the transport equation (I5.60p is solvable for e C"'^(0, T; "Hm) is due to 
the regularity of the eigenfunctions of the Stokes operators (cf. [3T] [IB])- Let 
_ ^™ ^ ^ C'^(0, T; "Hm) be the solution of the system of linear equations 

m J m m 

Y^AT\t)-gT{t) ^-Y.B2\t)hT{t)g^{t) -Y,ClgT{t)+V^-\t). 
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This system of linear equations is solvable because the eigenvalues of matrix of the 
coefficients ^1™"* (t) are bounded from below, since 

(5.61) AT^'y,-fj= / p\ijfdx>Mi \i;\^dx where V = TifZ-^- 

Jn Jn 

Thus we constructed a mapping M with A^(w™) — -u™. The energy estimate (|5.59p 
play the key to allow one to apply Leray-Schauder fixed point theorem for M. 

Remark 5.3. We know that from the estimates obtained in Section 2 and 3, we 
can pass to the limit and conclude that, there exists {p,u,d) such that, taking 
subsequence if necessary, 

p"" ^ p mLP{nx[0,T]) for anype (l,oo) 

M™ ^ M weakly in L'^{0, T; H^) and strongly in L2(0, T; L^) 

d"" weakly in L^{0, T; H^) and strongly in L^{0, T; H^). 

It follows easily from the above convergences that indeed (p, u,d) is a weak solution 
to the system (|l.ip . 



6. Appendix: Frozen Coefficient Method 

For completeness, in this section, we recall an application of the frozen coefficient 
method by considering the following problem (for more detail see |19j ) 

pvt — Aw + pv ■ Vw + Vp = /, 

V-t; = 0, i;(0)=0, v\9n = 0, 

in the domain 17 x (0, T). Here p is a given function in f2 x (0, T) satisfying 

(6.63) < Afi < p < M2 

(6.64) |Vp|<M3, \pt\<Mi 

(the derivatives of p may be the generalized ones). And / is a given function in 
L«(f7 X (0,r)) for q>^M■ 

The frozen coefficient method is now used to prove the following lemma. 

Lemma 6.1. If v & W^-^iiW^^iiVt)) and Vp G L?(17 x (0,r)) is a solution of 
problem h6.62]) . then for an appropriate constant a depending on q 

(6.65) \\vt\\L^ + ||Ai;||l, + ||Vp||l, < C(||/||l« + 11^^1120, 

where denotes the norm in L'?(i7x (0, T)), the constant C depends on q, Q, Mi, M2, M3 
and M4. 

Proof. Let (j^{x) be the set of nonnegative functions in C^(R") depending on the 
parameter A e (0, 1), forming a partition of unity in fl: 

k=l 

with 

suppC^ix) c 
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where the subdomains ft^ form a covering of fl which satisfy the foUowing properties 

(1) the diameters of $7^ do not exceed A; 

(2) there is a fixed number I such that the multiphcity of the covering of f2 by 
the subdomains Sl^ is < I; 

(3) Ef=i(Cfe'(a;))'>Ai>0,VxGf7; 

(4) \Di^(:^ix)\<cX-\^\ \p\ = i,2,yxen. 

For simplicity, we consider the case when p is independent of t. Fix Xk € f^^, for 
each k = 1,2, ...,N\, the vector Vk{x,t) = v{x,t)(^^{x) and the function pk{x,t) = 
p{x, t)(^ (x) are solutions in x (0, T) (in the same spaces as v and p) to the problem 

p{xk)vkt - A-Wfe + Vpk 

(6.66) = fCkix) + [p{xk) - p{x)]vkt - p{x)vk ■ Vwfc + 9k{x,t), 

V -Vk^v-VCkix), Vk{0)^0, Vk\on = 0, 

where 

n 

gk{x,t) = ~2Y,v.XL - vAC^+p\/Ci:+p{x)[vC{v ■ VC,^) + (C,^ - 1)C^ • V^;]. 

Assuming the right hand side of the first equation in (j6.66l) is the "free term" , 
making use of the results of Lemma 2.2 in [12], we have the estimate for Vk,Pk 

(6.67) \\vkt\\L-> + IIAwfclli. + llVpfellL. 

< C'ldl/lli, + max \p{xk) - p{x)\\\vkt\\Li + \\vk ■ Vwa;|1l<! + \\gk\\Li), 

where constant Ci depends on Af2. Since max \p(xk) — p(x)\ < M^\, and 

k 

the properties of C,^, we are able to derive estimate (|6.65l) from (|6.67p . by choosing 
appropriately A — (l+Ms)"^ min { 1, (2Ci)^^ } and applying the general Gagliardo- 
Nirenberg inequality (see ^). 

□ 

Remark 6.2. As seen in this simple example once we have the bound on the small 
oscillations the desired estimate follows. 
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